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Abstract 

Results of numerical simulation of the weak localization in two-dimensional systems in wide range 
of magnetic filed are presented. Three cases are analyzed: (i) the isotropic scattering and randomly 
distributed scatterers; (ii) the anisotropic scattering and randomly distributed scatterers; (iii) the 
isotropic scattering and the correlated distribution of the scatterers. It is shown that the behavior 
of the conductivity beyond the diffusion regime strongly depends on the scattering anisotropy and 
correlation in the scatterer distribution. 
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Interference (or weak localization) quantum correction to the conductivity arises from 
interference of electron waves scattered along closed trajectories in opposite directions. An 
external magnetic field applied perpendicular to the two-dimensional (2D) layer destroys 
the interference and suppresses the quantum correction. It results in anomalous negative 
magnetoresistance, which is experimentally observed in many disordered 2D systems.^ The- 
oretically, this problem was studied only for the case of random distribution of scattering 
centers and isotropic scattering (comparative analysis of different calculation approaches has 
been done in Ref. j^. As a rule, these conditions are not fulfilled in real semiconductor struc- 
tures. First of all, the scattering by ionized impurities dominates often in semiconductors at 
low temperatures. This scattering is strongly anisotropic, in particular, in heterostructure 
with remote doping layers. Besides, the impurity distribution is correlated to some extend 
due to Coulomb repulsion of the impurity ions at the temperatures of crystal growth. 

In the present work the role of scattering anisotropy in the quantum correction to the 
conductivity is investigated through a computer simulation. It is shown that the interference 
quantum correction is nonuniversal beyond the diffusion regime, when it is determined by a 
short closed paths with relatively small number of collisions. Its behavior with temperature 
and magnetic field change depends on the scattering details strongly. 

The weak- localization phenomenon can be described in the framework of quasiclassical 
approximation which is justified under the condition kpl ^ 1, where kp is the Fermi wave 
vector, / is the mean free path. In this case the conductivity correction is expressed through 
the classical quasipiohaJoiMtj density W for an electron to return to the area of the order 
Xpl around the start poinlj^*^^ 

5a = -ao — W = -2nl'^GoW, (1) 

TT 

where Go = e^/(2vr^^), = 27r/fci?, and ctq = nGokpl is the Drude conductivity. Prefix 
quasi- means that W takes into account the dephasing of interfering waves caused by ex- 
ternal magnetic field and inelastic processes. Expression allows ones to calculate the 
conductivity correction at any magnetic field and temperature.^ To find W and /, we simu- 
late the motion of a particle over the 2D plane with scattering centers in it. This technique 
was described in details in our paper, Ref. y. Here is outline only. The plane is represented 
as a lattice. The scatterers with a given cross-section are placed in a part of lattice sites with 
the use of a random number generator. A particle is launched from some random point, 
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FIG. 1: Sketch of model 2D system and fragment of the i-th. closed path near starting point. 

then it moves with a constant velocity along straight lines, which happen to be terminated 
by collisions with the scatterers. After collision it changes the motion direction. If the 
particle passes near the starting point at the distance less than d/2 (where is a prescribed 
value, which is small enough), the path is perceived as being closed. Its length and enclosed 
algebraic area are calculated and kept in memory. The particle walks over the plane until 
it escapes the lattice. As this happens one believes that the particle has left to infinity and 
will not return. A new start point is chosen and all is repeated. A fragment of one of the 
closed paths near the starting point is shown in Fig. ^ 

The simulation has been carried out for three systems. In system A the scatterers are 
distributed randomly, the scattering is isotropic that physically corresponds to a short-range 
scattering potential. In the system B the scatterers are distributed randomly and scattering 
is anisotropic. The angle dependence of scattering probability for this case is presented in 
the inset in Fig. Efb). It is close to that for the heterostructure with 6 doped barrier in 
which impurities are spaced from the 2D gas. The anisotropy shown can correspond, for 
instance, to a semiconductor heterostructure with impurity density of about 10^^ cm~^ and 
a spacer thickness of 50 A. Finally, the system C is characterized by isotropic scattering, 
but the distribution of scatterers is quite correlated. Namely, the nearest neighbor distance 
is not less than 20 whereas the mean distance between scatterers was about 30 (hereafter 
length and area are measured in units of lattice parameter and its square, respectively). 

The parameters used in simulation are the following: d = 5; the number of launches 
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FIG. 2: The interference quantum correction in zero magnetic field as a function of l/l^f,. (a) - 
Scattering is isotropic, the distribution of scatterers is random (system A). Dashed and sohd Unes 
are the analytical results, Eq. Q) and Eq. respectively, (b) - Scattering is anisotropic with the 
angle dependence of the scattering probability shown in the inset, the distribution of scatterers is 
random (system B). (c) - Scattering is isotropic, the distribution of scatterers is correlated (system 
C). (d) - The same data are presented in enlarged scale. Open and solid symbols correspond to 
Eqs. and (jSJ, respectively. 

10^.10^; the dimensions of the lattice is 10'' x lOl The cross-section of the scatterers is 
equal to 7. The density of scatterers is such that the transport mean-free path is about 120 
for all three systems. For illustration, let us set the lattice parameter equal to 5 A. Then our 
model provides an example of 2D system in which the transport mean-free path is / = 600 
A and Btr ~ 0.09 T. 
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As shown in Ref. |8|, the value of 5a for the model systems can be calculated as follows 



where summation runs over all closed trajectories among a total number of trajectories A^, 
b = B/B tr is the magnetic field measured in units of transport magnetic fields Btr = h/ {2eP), 
Si and li stand for the algebraic area and length of the i-th trajectory, respectively, is the 
phase breaking length connected in reality with the phase relaxation time through the 
Fermi velocity vp, = vpT^ (not to be confused with the diffusion phase breaking length 
Lfj, = y^Dr^, where D is the diffusion coefficient). Equation Q takes into account only 
the coherent backscattering correction to the conductivity. In order to take into account 
the nonbackscattering processes, which are important in the ballistic regime,-^ each term in 
Eq. Q should be multiplied by the factor [1 — cos(^i)],'^ where 6i is the angle between the 
first and last segments of the i-th closed path (see Fig. P): 



In what follows the results obtained from both Eq. ^ and Eq. will be considered. 

In Fig. 121 the results of our simulation for 6a{b = 0) are plotted against the ratio l/l^ 
(or that is just the same against t/t^). Open symbols correspond to the case when only 
the backscattering processes are taken into account, i.e., 5a is calculated from Eq. 
Data presented by solid symbols are obtained from Eq. (jHl) and, thus, include both back- 
and nonbackscattering contributions. In the assumption of oc T"^, where p > 0, this 
figure illustrates the temperature dependence of the interference quantum correction to the 
conductivity. The upper solid line in Fig. Efa) shows the known theoretical results for the 
backscattering contributions in weak localization^ 



The lower line is the results of Ref. |6| where both the backscattering and nonbackscattering 
contributions are taken into account 
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FIG. 3: The area distribution function of closed paths, obtained from the simulation procedure. 
Solid lines are the asymptotes for the purely ballistic and diffusive motion. Since W{S) is identical 
for S > and S < 0, only the range of positive algebraic area is shown. 

As is clearly seen the simulation and the theoretical data are almost the same. Thus, one 
can believe that the simulation procedure does not fail and for l/l^ > 3 x 10"'^ our model 
system is equivalent to an unbounded 2D system and the simulation gives correct results. 

Inspection of Fig. |2 shows that the temperature dependence of 6a in the systems with the 
anisotropic scattering (system B) and with the correlation (system C) is very close to that 
for the system with isotropic scattering and random distribution of scatterers (system A). In 
all three cases the temperature dependence of 6a is close to the logarithmic one for low l/l^ 
values, i.e., for low temperatures in real systems. The only difference is the absolute value of 
the interference correction for a fixed l/l^i, value. The relative difference in the temperature 
behavior for different systems becomes more pronounced at ///^ > 0.1 [see Fig. |^d)] when 
the quantum correction is mainly determined by the closed paths with a small number of 
collisions (so called ballistic regime). Thus, already analysis of the data carried out in the 
absence of magnetic field shows that the interaction correction feels the details of scattering. 
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Before to discuss the magnetoconductivity let us consider the area distribution function 
of closed paths W{S) because namely it determines the magnetic filed dependence of the 
interference induced magnetoresistance.™ As is seen from Fig. |21 the dependence 4:nl'^W{S) 
mostly follows S^^-lam for all three systems while S > l^. This is in agreement with the result 
of the diffusion theory for infinite 2D system which gives AttI^W{S) = S^^ tanh(7rS'//^) ~ 
for S > Thus, despite the fact that the distribution of scatterers is correlated in the 
system C, this system shows the diffusive behavior for large number of scattering events. A 
significant difference in W{S) is evident for small areas enclosed, S <C where the main 
contribution to W{S) comes from the closed paths with small number of collisions. It is 
obvious that details of single scattering event and the distribution of scatterers are important 
for statistics of short closed paths. For system A the area distribution function W{S) tends 
to the asymptotic S~^^'^ law valid when 5* -C [see Refs.QandQ. For other two systems, 
W{S) demonstrate more complicated behavior in our S'-range. 

So, this qualitative consideration shows that the low-field magnetoresistance {b <^ 1) 
which is mainly determined by the paths with large areas enclosed {S 3> must be similar 
for all three systems. Strong difference should appear only in high magnetic field (6^1) 
when short closed paths with S <^ P are suspended from the interference. 

Now we are in position to discuss the magnetoconductance Aa{b) = a{b) — o"(0). Namely 
this quantity is is usually measured to obtain the phase relaxation time experimentally. For 
this purpose an experimental Aa-vs-b plot are fitted to the one of the known theoretical 
expressions with as fitting parameter. The formula^ is widely used for this purpose 

Aa{b) = aGojv' + " ^ + ^) " ^^^j' (6) 
where 7 = t/t^, ip{x) is a digamma function, a is the prefactor. Theoretically, the prefactor 
is equal to unity. Experimentalist uses it as the fitting parameters together with (or 7). 
Below we follow this way considering our simulation data as the experimental ones. Then, 
we compare the fitting value of 7 with t/t^ used in the simulation procedure. 

Fig-EI^a) shows the low-field (6 < 0.2) magnetoconductance Aa{b) calculated from Eq. (jSJ 
for different t^/t values. As expected all systems demonstrate close behavior. It has been 
also found that taking into account nonbackscattering processes does not practically change 
the magnetic field dependence of A(j{b) for all the systems (insensibility of low- field magne- 



7 




FIG. 4: The low- field dependence of Act for different l^/l values (a) and the results of data pro- 
cessing of this curves with the help of Eq. © (b). 

toconductance to the nonbackscattering processes is discussed in Refs. O andO). 

If one considers the simulated Acr-vs-6 data shown in Fig. Ufa) as experimental ones and 
fits them by Eq. ^ as described above we obtain a nice coincidence in all the cases. How 
the fitting value of 7 matches the value of l/l^ = t/t^ put in Eq. is shown in Fig. Hl^b). 
As seen the fitting procedure gives the value of 7 which can differ by a factor of about two 
from the value of Z//,^. This may bee explained by the fact that Eq. (jH)) is derived within 
the framework of the diffusion approximation (7 <^ 1, 6 ^ 1), which is too strong.i^^ 

Figure El shows the dependencies of Sa{b), calculated for different t/t^ values with the 
help of Eqs. and in wider range of magnetic field including the ballistic regime 6^1. 
It is seen that the behavior of 6a {b) for b > 1 strongly depends on scattering details for both 
types of calculations. It is evident much better when considering the Acr-vs-6 data (Fig. Ej). 
So the magnetic field dependences of 6a and Aa are saturated in lower magnetic field when 
the scattering is anisotropic or distribution of scatterers is correlated. Such peculiarities 
directly follow from the difference in the area distribution functions at S < (see Fig. OJ. 
There is a appreciable deficit of the paths with small area enclosed, 5* < 200, in systems B 
and C as compared with system A. 

In conclusion, the magnetic filed behavior of the low-filed magnetoconductance in two- 
dimensional systems due to suppression of the quantum interference is not universal beyond 
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FIG. 5: The magnetic field dependence of 6a for different Z^// values, obtained with the help 
Eq. © (a) and Eq. © (b). 

the diffusion regime (6 ^ 1 or r/r^ ~ 1). It strongly depends on the scattering details. The 
use of theoretical expressions which were derived for the 2D systems with isotropic scattering 
and the random distribution of the scatterers for interpretation of the experimental data 
beyond the diffusion regime, that is a typical situation for experiments on high-mobility 2D 
heterostructures, can give inadequate information on the value and temperature dependence 
of the phase relaxation time. 
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FIG. 6: The magnetic field dependence of Act = (5(t(6) — (5(7(0) for different l^jl values, obtained 
using Eq. @ (a) and Eq. © (b). 
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